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Abstract 



We show how the standard domain wall action can be simply modified to 
allow arbitrarily exact chiral symmetry at finite fifth dimensional extent. We 
note that the method can be used for both quenched and dynamical cal- 
culations. We test the method using smooth and thermalized gauge field 
configurations. We also make comparisons of the performance (cost) of the 
domain wall operator for spectroscopy compared to other methods such as 
the overlap-Dirac operator and find both methods are comparable in cost. 
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I. INTRODUCTION 



Recently a great deal of theoretical progress has been made in the construction of lattice 
regularizations of fermions with good chiral properties [0-0] . For use in practical numerical 
simulations, though, approximations to these formulations are necessary. In the formulation 
using domain wall fermions (DWF) |I],|J, the extent of an auxiliary fifth dimension has to 
be kept finite in numerical simulations while chiral symmetry holds strictly only in the limit 
of infinite fifth dimension. The violations of chiral symmetry are expected to be suppressed 
exponentially in the extent of the fifth dimension [|]|| , but in practice the coefficient in the 
exponent can be quite small ||[7J and the suppression correspondingly slow. 

In the case of overlap fermions there is no such problem in principle. However, there is 
a problem of practicality: how to deal efficiently with e(H) = H/ y/ H 2 , where H is some 
auxiliary Hermitian lattice Dirac operator for large negative mass, but free of doublers. Most 
commonly, the Hermitian Wilson-Dirac operator H w is used. So far, the best methods use 
rational polynomial approximations of t(H) rewritten as a sum over poles 

Variants of this are Neuberger's polar decomposition || and the optimal rational polynomial 
approximation of Ref. j^JlOH. In the application of e(H) on a vector ip the shifted inversions 
(H 2 + bk)~ lr 4) are done simultaneously with a multi-shift CG inverter [IT}]. This is referred 



to as the "inner CG", since typically another, "outer", CG method is used to compute 
propagators or eigenvalues of the overlap Dirac operator containing e(H). 

In these methods, the main problem is the difficulty of getting an accurate approximation 
to e(x) for small afj. As emphasized in Ref. flT2"|JI0|| accuracy for the lower range of the 
spectrum of H can always be enforced by calculating a sufficient number of eigenvectors of 
H with small eigenvalues and computing e(H) in the space spanned by these eigenvectors 
exactly. The approximation is then only used after projecting out these eigenvectors. This 
is referred to as "projection". 

Small eigenvalues of H w occur quite frequently on (quenched) lattices used in current 
simulations. Indeed, a non-zero density of eigenvalues near zero has been found jl3|| . Pro- 
jection is thus essential in efficient implementations of the overlap Dirac operator. It is 
known that the zero eigenvalues of H w correspond to unit eigenvalues of the transfer matrix 
T that describes propagation along the fifth direction for domain wall fermions 0. The 
left and right handed physical fermions in the domain wall formulation occur along the two 
boundaries in the fifth direction. A unit eigenvalue of the transfer matrix then allows for 
unsuppressed interaction between left and right handed fermions and thus to a breaking 



1 Using the scaling invariance e(sH) = e(H) for any positive scale factor s the upper range of the 
spectrum of H can always be scaled to be in the range of good accuracy of the approximation to 
e(x). 
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of chiral symmetry. Therefore the same modes that make the approximation of e(H) diffi- 
cult to achieve in overlap fermions, and need to be projected, are responsible for the chiral 
symmetry breaking for domain wall fermions at finite L s . 

If we could project out the modes with near unit eigenvalue of the transfer matrix and 
take their contribution with effectively infinite L s we could obtain a domain wall fermion 
action with no chiral symmetry breaking even at finite L s . This would make domain wall 
fermions with finite L s equivalent to overlap fermions. It would then become a matter of 
computational efficiency as to which form is to be preferred for numerical simulations. Based 
on ideas of Borigi |TJ[ we show in section |TJ how projection for domain wall fermions can be 
achieved. In section [TIT] we illustrate how projection works, first on smooth instanton con- 
figurations and then on a configuration from a quenched simulation. We discuss the degree 
of chiral symmetry conservation that can be achieved and compare the performance (cost) 
of domain wall fermion implementations with and without projection to overlap fermions. 
We conclude the paper with a brief summary and some discussions in section |V|. 

II. PROJECTED DOMAIN WALL FERMION ACTIONS 

Kaplan's method M realizes a single massless fermion field through an infinite tower (or 
infinite number of flavors) of massive fermion fields with a particular flavor structure. The 
flavor index can be interpreted as an extra (here fifth) dimension, and the flavor structure 
as a defect along the fifth dimension to which the massless fermions are bound. Lattice 
calculations necessarily require a finite fifth dimensional extent and lattice spacing, as well 
as regularizations for the derivatives in the four dimensional space and in the fifth direction. 
These regularizations are not unique. In particular, higher order terms in the fifth dimen- 
sional lattice spacing can be added to the action, as long as they do not change the "defect 
structure". Such changes can affect the effective four dimensional theory by altering the 
discretization effects in the four dimensional lattice spacing, but not the continuum limit. 

The finite fifth dimensional extent induces some chiral symmetry violation. This can be 
eliminated by "projecting" an appropriate number of violation inducing states. We show 
how the projection method can be implemented for two different five dimensional domain 
wall actions. The Wilson fermion action is used for the four dimensional part of the action, 
but it should be understood that other actions realizing a single massive fermion field could 
be considered, and much of the following derivations do not depend on the specific form 
chosen. 

A. Form of the actions 

To fix our notation/conventions, we write the usual, 4-d, Wilson-Dirac operator as 

, 4 

D W (M) = (4 + M)S X , V - - [(! " J»)U,{x)5 x+ ,, y + (1 + 7„)^(</)<W 
_(B + M C \ 
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Here, with = (ak,il), 

- 4 

2 [ 2<J *.» ~ c/ m( x ) - u l(y) s »,v+n ■ ( 3 ) 

We are using a chiral basis: 7 M = , 75 = . We will also use the hermitian 

Wilson-Dirac operator 

H W (M) = l5 D w {M) = + M _ B °_ M ) . (4) 

Usually, we will use a large negative mass here, but often omit the argument from D w and 
H w . 

In this notation, the usual domain wall fermion action of Shamir || reads, with an 
explicit fifth dimensional lattice spacing 05 for the hopping term — the 4-d lattice spacing 
a is kept fixed at a = 1 throughout — , 

l s 

S DW = -VDf w i2 = -Y d %{ [a 5 D w (~M) +1]%- P^ i+1 - P+%-i] (5) 



i=l 



where the extent of the fifth dimension, L s has to be taken as even. The gauge fields in 
D w (— M) are independent of the fifth coordinate, and the fermion fields satisfy the boundary 
conditions in the fifth direction: 

P-*l.+i = -mP-Vi, P + ^ = -mP + ^ Ls . (6) 

P± are the chiral projectors, P± = |(1 ± 75). < m < 1 is proportional to the quark mass 
(for small m). The 4-d physical fermion degrees of freedom are identified with the fields at 
the boundaries as 

q R = P + V Ls = q L = P_^i = ttf , 

q R = y La P_ = yl, q L = ^P + = ^. (7) 

As will be shown below, projection of low-lying eigenvalues of — logT, with T the transfer 
matrix along the fifth direction, can be achieved by introducing an additional term in so 
that the complete domain wall fermion action reads 

Sdwp = -*£>gk P * = ~ M-(~ M ) + ~ " 

i=i 

+ M(m)[P_*i+P+* L J. (8) 

Borigi |14]] introduced an interesting variation of the domain wall fermion action, which 
differs from the usual one by terms of order 0(0,5). He called it (5-d) truncated overlap 
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action, since its effective 4-d action for the light fermions is just the polar decomposition 
approximation of order L s /2 to the overlap Dirac operator introduced by Neuberger ||. In 
our notation, and introducing again an additional term to be used for projecting low-lying 
eigenvalues, Borigi's variant reads 



S DW , = -VD® W ,V = -^2%{[a 5 D w (-M) + l]^i+[a 5 D w (-M) - l]P-V i+1 + 

i=i 

[a 5 D w (-M) - + f !i(m) [P_*x + P+V Ls ] ■ (9) 



The boundary conditions in the fifth direction remain as in 
operator is, including the boundary conditions, 



The kernel of the 5-d 



D 



(5) 
DW 



AP_ D_P_ 



D-P + 




D-P+ 








P_P_ 
D, D_P_ 

























D + 
DP 



-mD_P + - 





apA 



(10) 



where 



D ± = a 5 D w (-M) ± 1 



(11) 



The kernel of the standard domain wall operator, but including the term to be used for 
projection, D^^p, is simply obtained by the replacement D_ — > — 1. 

We now integrate out the heavy fermion degrees of freedom to arrive at a 4-d Dirac 
operator describing the light fermions. Following Ref. |14| we define V by 

, L s - 1 



PS jjk + P+S j+1>k for j = 1, 
P-^L a ,k + P+Si, k for j = L s 



(12) 



This has an inverse V 1 = V', given by 



i jk 



P_6 jjk + P + 6j_ 1>k for j = 2, 
P_5i, fc + P+5 La:k for j = 1 



(13) 



Next, we introduce Xi s through = (Vx)i an d define 4-d operators Q± as 



a^H w P± ± 1 for the standard domain wall action 
asif,,, ± 1 for Borigi's domain wall action . 



(14) 



Then, using 7sP + = P + and 75P- = — P_ as well as the boundary conditions on the fermion 
fields, we can rewrite both domain wall actions as 
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S (5) = - {*i>y s [Q-(P- - mP + ) X i - idxi + Q+X2] (15) 

Ls-l 

+ J^^MQ-Xi + Q+Xi+i] +^lMQ-Xl s +Q+{P+-mP_)xi]}. 

i=2 

Next, we introduce = XiQ- lr Y5 an d 

t- 1 = -g: 1 g+ . (16) 

This change of variables would give rise to a Jacobian in a dynamical simulation. However, 
it would be exactly canceled by the Jacobian of this transformation for the pseudo-fermion 
fields which are needed to cancel the bulk contribution from the 5-d fermions. The action 
for the pseudo-fermions is obtained by the replacement m — > 1 from the fermion action. The 
Jacobians cancel since integration of fermions (Grassman fields) acts like differentiation. We 



note that for the standard domain wall action Q± do not commute and the ordering in (16) 
is important. From (|J) we find 

1 a 5 C±\ / B 0\ (B + alC^ a 5 Ci 



where B = 1 + a^,(B — M), and thus 

T is the usual domain wall fermion transfer matrix || in our conventions. 

For later use it will be convenient to introduce a 4-d Hamiltonian Ht @] such that 



= ' (19) 

For Borigi's domain wall action we have simply Ht = H w , while for the standard domain 
wall action one finds 

Ht = H w = H w . (20) 

2 + a 5 H w ^ 5 2 + agysHw 

In terms of the new fields x an d X t ne 5-d actions become 
= - X Dfx = ~ I Xi [(P- ~ mP+) X i - Q--lsA X i - T~\ 2 ] 

+ Yl * bd - T ~ 1 X l+ i] + Xl 3 [xl s - T~\P + - mP-) X i] V (21) 

i=2 ' 



Now, we integrate out, in succession, Xl s ,Xl s , Xl s -i, Xl s -i, ■ ■ ■ , X2, X2- For this we use, 
at the first step, 
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XuXu ~ Xl s -iT l x La - XlT \P + - mP_) X i = 

[xl s - Xl s -iT- 1 ] [ XLs - T-\P + - mP-) X i] ~ Xl s -iT~ 2 (P + - mP_) Xl 

and at the (L s — i)-th step 

Xi+iXi+i ~ X^Xi+i - X,+iT- Ls+ \P+ - mP-) X i = (22) 
[Xi+i ~ ^T- 1 } [Xi+i ~ T- L ° + \P + - mP_)xi] - Xl T- Ls+l -\P + - mP-) X i- 

With a change of variables Xi+i = Xi+i ~ T~ Ls+t (P + — mP_)xi and x'i+i — Xi+i ~ XiT' 1 
the integration over Xi+i->Xi+i is trivial, giving a factor 1. At the end, we arrive at the 4-d 
action for X i, X i 

S {4) = -Xi [(P- ~ mP + ) - T~ L °{P + - mP_) - QzH^Xi = -Xi^ (4) (^)Xi- (23) 
The kernel is 



D {i \m) = (P_ - mP + ) - T- Ls (P + - mP_) - QZ l ^A 

—— (T" L * + 1) 75 + —z- {T- L ° - 1) + QZ l i,A 



(24) 



{T-Ls + X ) 7s 



1 + m 1 — m T 



+ 



"75 



T-Ls + 



Finally, integrating out X ii Xi an d dividing by the pseudo-fermion determinant, obtained 
by the substitution m — > 1 and requiring A(m = 1) = 0, we find 



det£>( 5 )( 



m) 



^w=d*{[ D «(i)]- , ^H} 



detP( 5 )(l) det. 
Now, from eq. (f?|), D (4) (l) = - (T" Ls + l) 75 and thus we get 



(25) 



[D^\l)Y l D^\m) = - 2 

Using flT9| ) we note that 
t -l s _ 1 



1 + m + (1 — m)7 5 



T + 27 5 [Q_(T- L * + 1)] SsiH 



• (26) 



;i + a 5 if T ) L s-(l-a 5 # r ) L * 



T- L * + 1 (1 + a 5 # T ) Ls + (1 - a 5 H T ) L ° 



(27) 



where £ n (x) is Neuberger's polar decomposition approximation to e(x). Ignoring the term 
with A(m) for the moment, we see that 



1 + m + (1 - m)j 5 e Ls / 2 {a 5 H T 



D to Jm) . 



(28) 



This is just Neuberger's polar decomposition approximation || to the overlap Dirac operator 
for auxiliary Hamiltonian H?, which we denote by D tov , "truncated overlap." For L s — > 00 
it becomes the exact overlap Dirac operator D ov Alternatively, we can write 
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T- 



tanh ( — — log |T 



(29) 



and eq. (|26"D, still ignoring the term with A(m), becomes the effective 4-d Dirac operator that 



Neuberger derived for domain wall fermions in [16[ . Here we have used that L s is even and 
written the formula in terms of the absolute value |T| to indicate that everything remains 
well defined when an eigenvalue of T becomes negative. 

We now can use A(m) to project out low-lying eigenvectors, Wj, of the auxiliary Hamil- 
tonian H T for which e La / 2 with finite L s is not a sufficiently accurate approximation to e(x). 
Let 



H T Vi = XiVi 



where, from eq. ([19]), 



Tvi = TiVi, 



1 - a 5 \i 



Pi = ViV] , 



1 + a 5 \i 

The projection can then be achieved by setting 

A(m) = (1 - m)7 5 Q_y^g i P i; 



(30) 



(31) 



(32) 



with 



9i = 2 [- ( T r Ls - !) + ( T r Ls + 1) e^Ai)] 



(33) 



Note that A(m) vanishes for m = 1, as required for the pseudo-fermions. With this A(m) 
we find, instead of (BSh, 



[,D( 4 )(1)] 1 J D( 4 )(m) = J D OT (m) = 

-| 1 + m + (1 - m) 75 \e Ls/2 (a 5 H T ) < a ^)^ 



(34) 



For Borigi's variant of the domain wall action, where Ht = H w commutes with Q± the 
"projection operator" A(m) can be simplified to 



i(m) = (l-m) 75 ^/ i P i; 



(35) 



with 



ft = - 1) [- (Tf Ls - 1) + (Tr Ls + 1) e(a 5 A„ 

J(a 5 A<-l) for Ai > 



(36) 
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B. Propagator 



Our next step is to relate the 4-d overlap propagator to the 5-d propagator of the corre- 
sponding domain wall fermion action. In all steps below A(m) for the eigenvalue projections 
is included. From fl34l ) we find, obviously, 

D- v \m)=[D^(m)]- 1 D^(l). (37) 

To connect this to the 5-d theory, we consider, motivated by the fact that the light 4-d 
fermion is q = (V^ 1 ^)! = xi> 

X = jp- 1 [D( 5 >(m)] -1 D (5) (1)P} 



-5(5) 



(38) 



7^7 / \{dXidXiY.^Q- 1 ^ [^CO^L* 

J i k 

^/n rf X^£x^?(l)] fcl e- S(,> . 



5(5) 



Here we have used the 5-d = Qz 1 "f5D^V introduced in eq. (pT|). Now, we integrate 
out, in succession, XL s ,X,L a , Xl s -i,Xl s -i, ■ ■ ■ , XiiXi- From the transformations ( p2[ ) we see 
that Xk XiT~ k+1 in the process. Therefore, we obtain 



»0 (4) (m)xi 



But from pip we see that 



ki 



[P--P+ 



T- L °(P,-P_) 



(39) 



(40) 



where we used the fact that A(l) = 0. Thus we finally obtain 



)(4), 



in 



1 D^\l)=D- ov \m) 



To solve 

D ov (m)tfj = b 
we introduce b = (6, 0, ... , 0) T and solve flip 

D (5) (m)(j) = D {5) (l)Pb. 

ip is then obtained as 

1 p=(r- 1 <p) 1 . 

The physical propagator has a contact term subtracted |12|] , and we arrive at the general 
form for both auxiliary Hamiltonians considered here, with and without projection, and also 



(41) 

(42) 
(43) 
(44) 



valid for finite L, 



D- V \m) 



1 



1 — m 



[D-\m) - 1] = \V' X [£> (5) M] 1 D®(1)V\ - 1 

L J 1 — m LI L J J li J 



(45) 
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C. Relation of the 5-d and 4-d operators 



So far, we have related D ^(m) and det D ov (m), or their truncated versions, to 5-d 
operators, see ([25]), fl28|) and (0) and finally pl|). E.g. for eigenvalue calculations it 
would be useful to establish a similar connection for D ov (m) or D tov (m). Consider D^\m) 
introduced in (|2ll). We find 



[Df{m)- 
1 - m)P, 



4 5) (1)^ 
- QZ 1 j 5 A( m ) 



(46) 



5 il 5 jl -T- 1 (l-m)P^S iLs S 



Thus 



[4 5) (1)] Df{m) = 

1 + l D ?0)]a l^ 1 ~ m ) P + " Q-^Mm)} S fl - (1 - m) [Df{l)]- l La T-'PS^ 



(47) 



So only the first column (of 4-d blocks) is different from the 5-d unit matrix. Let's call the 
entries in the first column X±, X 2 , . . . , Xi . Similarly, its inverse is 



i - [nm ( 



b. (5) (i; 



(48) 



il 



'1 - m)P+ 



5n + (1 



m 



[Df{m 



iL, 



T- l P-5 3l 



Again, only the first column is different from the 5-d unit matrix. Let's call the entries in 
the first column Y\, Y2, . . . , Yl s - Since these 5-d matrices are the inverse of each other, their 
product is the 5-d unit matrix, i.e. 



Y 1 X 1 
Y2X1 + X2 



1=>Y 1 = (X,)- 1 
0^Y 2 = -X 2 (X 1 ] 



YjXi +Xj = Q=> Yj = -X j (X lj 



for j 



The first equation, Y\ = (Xi) 1 establishes the relation we are looking for 

D ov ( m ) = [ J D( 5 )(l)] _1 J D (5) (^)^ p } ii • 



(49) 
(50) 



(51) 



Note that the components of the 5-d propagator used for the inverse of the overlap Dirac 
operator in eq. (41) do not correspond to the physical quark propagator (qq) as obtained from 
domain wall fermions. The reason is that the overlap propagator still needs a subtraction 
of a contact term and a multiplicative normalization jll|. We consider here the standard 
domain wall action without projection, and denote the corresponding 4-d Dirac operator 
Dtov(™>', Ht) ■ To make the connection between the subtracted 4-d propagator and (qq) 
explicit, we write the 1 for the subtraction as 



1 



v- 1 



D%(m) 



-1 



D%{m)V 



(52) 
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and thus 



D^(m:H T )-l = \v- 1 



-i 



)(5) ( 

'dw\ 



)(5) 

dw y 



V 



(53) 
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From the matrix for D^ w , eq. (|T0|) with the replacement D_ — > —1 and setting A(m) = 0, 
we see that 



)(5) ( 



)(5) 



1 - m) [P-^^i + P + 5n5 jLs 



Therefore we obtain 

D~l(m;H T ) - 1 = (1 - m) {rl^H]- 1 ^}^ 



(54) 



(55) 



where = <5i,L s +i-j is the inversion operator of the fifth direction. Now from (^) we see 
that the physical fermion degrees are given in terms of the domain wall boundary fermions 

as 



Hence we find 



[P~ 1 [D^ w (m)r 1 JP} n = ^ [D^(m;H T ) l] 



(56) 



(57) 



The usual domain wall physical fermion propagator automatically contains both the sub- 
traction and multiplicative normalization of the overlap fermion propagator. 

The hermitian conjugate of the 5-d operator Dp W ,, eq. (|l0l), of Borigi's variant is different 
than that of the usual DWF operator, D DW) without projection. In particular it does not 



have the (generalized) 75 hermiticity, D^^ = D K ^ W ^^J , with J the inversion operator 
of the fifth direction. Instead we find 



)(5) 



D 



(5)t 
DW 



\-mP + Dl - P+At 



P+DL 
P_Dl_ 










p+dI 

Di P+Dl 















-mP_Dl\ 















P_D ] _ D\ 
P_D ] _ 



P + D ] 



(58) 



where 



Dl = a 5 Dl(-M) ± 1, A\m) = (1 - m) ^ 



7s 



(59) 
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D. Eigenvalues of the pseudo-fermion operators 



Here, we want to investigate a little more the properties of the pseudo-fermions, needed 
to cancel the bulk contributions in the 5-d domain wall fermion approaches. 

For both the standard domain wall fermion action @ and for Borigi's variant (|]) we 
find for the pseudo-fermion matrix, including the boundary conditions and recalling that 
the projection operator A{1) vanishes, 



Z>(6)(l)Z>(5)t(l) 



( X Y •■ • 
Y X Y ••• 
Y X Y ••• 



-Y\ 






-Y 



Y X Y 
Y X 



(60) 



with X and Y 4-d hermitian matrices, 



Y = -P + D, 
X = D+D\ + 1 

for the standard domain wall action, and 



(a 5 D w + a 5 Dl + 2) 



2Y 



(61) 



Y 
X 



D+P+Dl 



D_P_D\ 



alD m Dl - 1 



2a\D w Dt + 2 = AalD w Dl 



2Y 



(62) 



for Borigi's variant. Similar relations can be found in Ref. [17 



Now, let S be the shift (or translation) operator in the 5-th direction with anti-periodic 
boundary condition: 



/ 1 
1 



S 



o o\ 






-10 



1 
0/ 



(63) 



It is easy to see that [^^^(^^^(l)] = 0, and so S and D® (l)D^ (1) can be diaeo- 
nalized simultaneously. The eigenvalues of S are exp{in(2k + 1)/L S } for k = 0, . . . , L s — 1, 
because of the anti-periodic boundary conditions. The corresponding eigenvectors are 



w(kf = 

(t> 4 , t> 4 exp{in(2k + 1)/L S }, t> 4 exp{z7r2(2fc + 1)/L S }, . . . , v 4 exp{i7i(L s 



(64) 

l)(2k + l)/L s }) T 



with t> 4 some 4-d vector. This also has to be an eig envector of D (5) (l)D (5)t (l). From 
and using flBT]) and (|62| ) we obtain the eigenvalue equation for w 4 
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Ja 2 5 D w Dl - 2Y 



cos I ^-{2k + 1) 



V4 = \(k)v4 



(65) 



Here J = 1 for domain wall fermions, and J = 4 for Borigi's variant, and we have indicated 
the dependence of A on k, the momentum in the 5-th direction. Since the left-hand-side 
of (|65|) remains unchanged under k — > L s — fc — 1 we conclude that A(L S — A; — 1) = X(k). 
Therefore, the eigenvalues of D^(1)D^(1) are all (at least) doubly degenerate. 

Furthermore, for large L s and small k the pseudo-fermion domain wall eigen- 
values closely track the Wilson eigenvalues, the eigenvalues of D W D\ J01 since then 



2Y 



1 - cos ( -H2A; + 1) 



rs Yir 2 (2k + l) 2 /£ 2 m (H) is only a small perturbation. In 
particular, at zero crossings of the Wilson Dirac operator, a pair of pseudo-fermion domain 
wall eigenvalues will go to zero for L s — > 00. 



From the domain wall Dirac operator Z? 



can make a hermitian version, 



(5) 
DW 



[m) without the projection matrix A(m), we 

the inversion 



^SJy( m )^75, with J^- 



operator of the fifth direction. From eq. fllUD with P_ — > — 1 and A = it becomes 




















P_ D+75 \ 
P_ D +75 -P + 
P +7 5 -P+ 



P_ P +75 -P + 
\P +75 -P+ 















— mP_ 



(66) 



are 



Since P^]y(m)P^|^(m) = [ZZ{^(m)] 2 , the vectors u>(&) and w(L s — k — 1) of eq. 
degenerate eigenvectors of [P^^(l)] 2 . To get the eigenvalues of ZZ^^l) we need to compute 
the 2 x 2 matrix (w(k 1 )\H^ ) w (l)\w(k 2 )) with k u k 2 = k, L s - k - 1. We find 



<MM/41(i)H^)> 



g e _ ;& (2 fcl+ i)0--i) e ^(2 fc2+ i)(L.- J -) (u4 | P _ e -^(2* a +i) + D+J5 _ p+e^+i)^ (67) 

The sum over j gives the constraint k 2 — L s — k\ — 1, and we thus find the 2x2 matrix to 
be of the form 

b 
b* 

and therefore the eigenvalues of the hermitian pseudo-fermion domain wall operator come 
in ± pairs. 



E. Preconditioning 



Even-odd preconditioning is a common technique to reduce the condition number of a 
matrix |T3] and can be implemented in both versions of the domain wall actions discussed 
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above. With projection, the matrices have a more complicated form reminiscent of precon- 
ditioning for the clover fermion action. 

Write the matrix as a two by two block matrix 



Aee Be 



EO 

oo 



(68) 



where the upper case characters E and O label even and odd (checkerboard) sites, respec- 
tively, in five dimensions. The matrix can be brought into an even-odd block diagonal 
form using the lower and upper block triangular matrices 



l EE 







BqeA ee Iqo 



U 



Aee Beo 
loo 



with the transformation 



m = l-'d^u- 1 



Iee 








Aoo — B eA ee B E o 



(69) 



(70) 



This decomposition always exists as long as Aee is non-singular. 

The structure of Aee is deduced from the structure of D W (M) and the projection matrix 

(5) " (5) 

A that occur in both the standard domain wall and Borici variants, D DW and D DW ,, respec- 
tively. In Borigi's variant the nearest neighbor coupling in the fifth direction contains also 
hoppings in the spatial directions leading to a complicated even-odd structure. This case 
will not be discussed further here. For preconditioning of the standard domain wall action 
with projection, after a rescaling of the fermion fields to bring D^ w into the "kappa" -form, 
i.e., having 1 in the diagonal except for the projection matrix A, we have 



/ 1 - i ee P_ 



A 



EE 



■ -A eo P + \ 

■ 



V 













• 1 

• 1 



(71) 



/ 



Here A is rescaled by 2k, with k = 1/(2 (05 (4 + M) + 1)), and e and o refer to the even and 
odd sub-lattice in the 4-d checker-boarding. 
For even-odd preconditioning, we need 



f(i-A ee p^y l 
1 



A' 1 



EE 



(I- A ee P^)- l A eo P + \ 




V 













1 







1 



(72) 



/ 



Therefore, for every application of the projected domain wall Dirac operator, the inversion 
(1 — A ee P^)~ l (j) for some 4-d even sub-lattice fermion has to be performed once — the 
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two occurrences in A EE can be combined into the application on one 4-d even sub-lattice 
fermion field. In the limit of the fermion mass m — ► 1, the matrices A EE and Aoo become 
the unit matrix. For small m, the matrix (1 — A ee PJ) should be well conditioned, and the 
cost of inversion is the number of iterations involving applications of D w on half of the four 
volume. This "inner" iteration cost is a multiplicative overhead on the cost of inverting D^ 5 \ 



F. Dynamical fermions 

Implementing an HMC algorithm for dynamical domain wall fermion simulations is a 
straightforward modification of dynamical Wilson fermion simulation code. Implementing 
projection does not cause any fundamental difficulties. What is needed, in addition to the 
usual case, for the computation of the fermion contribution to the force is dA(m; H(U))/dU, 
where we indicated the implicit dependence on U of the projection operator A. Using the 
chain rule, what we need is 

du and m = m vl+Vi m> (73) 



for H(U)vi = \{Vi with v\vi = 1. We easily find 

d\i *dH(U) 



dU 1 dU 



-Vi ■ (74) 



and 



dU v ; 



' f dH(U) dHOJ) 

tl- —!)■ — — 

1 dU 1 dU 



(75) 



Therefore, a different 4-d inversion is needed for every projected eigenvector i^. For the 
standard domain wall fermion action, where H = Ht, additional 4-d inversions are needed 
to compute ■ 

We note that the modifications involving projection for the force term may not be neces- 
sary in practice. The unprojected action can be used to construct the guiding Hamiltonian 
for HMC with the projected action used for the computation of the initial and final energies 
in the Metropolis accept/reject step. A critical factor in the choice will be the acceptance 
rate which will suffer if the guiding Hamiltonian is not accurate enough without projection. 



III. RESULTS 

A. Smooth gauge field 

We start our analysis by first considering the case of spectral flow on a smooth SU(2) 
instanton configuration. We use the instanton construction from ref. |TI|, namely a single 



3 4 instanton with p = 1.5 and anti-periodic boundary conditions. In this section we will 
mainly focus on the standard domain wall operator. Similar results are found for the Borigi 
variant. 
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What we are interested in is first determining how various domain-wall like actions 
reproduce topology on the lattice. Without projection, we note that the index of the massless 
four dimensional (4-d) operator is 

Q = Tr( l5 D tov ) = ^Tr (e L , /2 (a B H{-M))) (76) 

with the auxiliary Hamiltonian H given by either Ht or H w — the two cases under consider- 
ation. For infinite L s , El s /2 becomes e and Q is a measure of the discrepancy of positive and 
negative states of H. A simple way to determine Q is to start from some M where we know 
that Q = and increase it to the desired value while counting the change of the number 
of positive states. We emphasize that this procedure works as long as the spectral flow is 
smooth — namely, as long as the eigenvalues of H are well behaved. While H w is always 
well behaved, we note that this is not the case for the domain wall Ht(—M) for M > 2. In 
particular, the spectral flow method will fail for M in the "doubler" region where there are 
multiple species of chiral fermions. We will address this point in more detail later. 

In Fig. [l| we show the spectral flow of the lowest 10 eigenvalues of the hermitian Wilson- 
Dirac operator H w (—M), the lowest 10 eigenvalues of the hermitian domain wall Dirac 
operator J' , y 5 D^ w (0) with a 5 = 1 for various fifth dimensional extent values L s , and the 10 
lowest eigenvalues of the hermitian overlap-Dirac operator (with projection, i.e. with chiral 
symmetry), all as function of the "domain wall height M". We recall from eq. ([51]) 

D ov (0) = {V- 1 [D^{l)]~ X D^{{S)v) n (77) 

that when D ov (0) has zero-eigenvalues we can expect D^ 5 \0) to have zero-eigenvalues as 
long as D^(l) is not singular. We see in Fig. [I] that there are exact zeros of D ov (Q) for 
M beyond the crossing, and the spectrum changes discontinuously at the crossing |§. For 
finite L s , the D^^(O) eigenvalues decrease quickly in L s . For M beyond the crossing they 
quickly go to zero and in fact appear exponential in L s , but near the crossing, when the 
eigenvalues of H w are small, the decrease is slowed. 

This slow rate of convergence in L s , when the eigenvalues of H w are small, is a prime 
source of difficulty in recent numerical calculations involving domain wall fermions. The 
reason for the slow rate of convergence is clear when the form of the 4-d operator is examined, 



D t0 v{m) = - 



l + m + (l-m)^e Ls/2 (a 5 H) . (78) 



When SL a /2{a^,H) deviates from one, for example for the smallest and largest eigenvalues 
of H, there is violation of chiral symmetry. Fig. |2] shows the domain wall spectral flow 
when the 5 smallest (closest to zero) eigenvectors of Ht are projected out as prescribed by 
eqs. (p2|),(^3|) and (|34l) . With projection turned on, D^^(0) no longer has (generalized) 75 



hermiticity, so the eigenvalues of +y (D^ w (0))W^ W (0) are determined. We see that just 
beyond the crossing there is a discontinuous drop of the eigenvalues to zero. In fact fairly 
small eigenvalues are seen even for L s = 4. Also shown are the lowest eigenvalues of Ht- 
We see they agree quite well with the eigenvalues of H w up to a factor of two as predicted 
byeq. (EO). 
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All computations of eigenvalues and eigenvectors where done with the Ritz method J20| . 



For computing eigenvectors of Hj<(—M), in each Ritz iteration, the application of 
Ht{—M)^Ht{—M) on a vector is needed. From eq. ( p0[ ) for Ht, the expression for H T H T 
can be combined into 

H T H T = l5 D w - - l5 D w , (79) 

(2 + a 5 Dl)(2 + a 5 D w ) 

so a single inversion of a hermitian positive definite operator is needed for each Ritz itera- 
tion and Conjugate Gradient for Normal Equations (CGNE) was used. Typically 30 to 60 
iterations of CGNE were needed for each inversion to achieve an accuracy of 10 - '. 

It should come as no surprise that zero eigenvalues of the domain wall operator can be 
obtained when projection is turned on since this is the same mechanism by which they are 
obtained for the overlap-Dirac operator |§. What was key to the latter case was enough 
projection coupled with a judicious choice of the fifth dimension spacing 05. Since a 5 only 
serves as a multiplicative scaling of H w , as can be chosen so that the lowest unprojected 
and highest eigenvalues of H w lie within the range of approximation of e. For the overlap- 
Dirac operator in Ref. 0, an optimal rational approximation was chosen; however, this 
merely increased the useful range of the approximation. When the condition is met for a 
valid approximation of e, there are no other cutoff effects associated with having a finite 
05. For example, observables of the effective 4-d theory will not have slightly different 4-d 
lattice spacing dependencies for slightly different choices of 05. However, for the conventional 
domain wall action, Ht intrinsically depends on 05 as seen in eq. (0). Therefore, even for 
infinite extent in the fifth direction, where e(a$HT) converges to e(asifr) = c(Ht), choosing 
different a 5 results in different 4-d lattice spacing dependence for physical observables. This 
fact limits the usefulness of adjusting a 5 with the goal of imposing arbitrarily precise chiral 
symmetry for (small) finite fifth dimension as will be shown later. 

In Fig. |3] we examine more closely the smallest eigenvalue of the domain wall operator 
and the deviation of e(a 5 Hj<) from one for a§ = 1. We see that without projection the 
eigenvalues decrease approximately exponentially in L s with a slow rate varying with M, 
while with projection they drop to about 10~ 4 after the crossing. This is the accuracy to 
which the eigenvalues were computed in single precision and is "zero" here. For increasing 
M, the zero eigenvalues slowly increase. 

Some clue to the slow increase can be seen from the right panel in Fig. |3| which shows 
the deviation from unity of e(Aj) where Aj is the first, fifth and largest eigenvalue of Ht- 
The deviation for the fifth and the largest eigenvalue gives a measure of the accuracy of 
the approximation since the first five eigenvalues are projected out of Ht- The deviation of 
the first eigenvalue indicates how much the standard unprojected case deviates. For L s = 4 
the deviation of the largest eigenvalue is clearly visible, and for L s > 8 the deviation is 
off the bottom of the graph until M approaches 2. For L s = 8 a typical deviation of the 
fifth eigenvalue is about 1% while nonetheless a domain wall eigenvalue of 10~ 4 is obtained. 
However, as M is increased, the largest eigenvalue is not well approximated and the domain 
wall eigenvalue deviates away from zero due to the resulting chiral symmetry breaking. 

We will conclude this portion of the analysis with a final look at the role of the pseudo- 
fermion term in eq. (|51|) (repeated again in eq. (|77|)). It was stated that zero eigenvalues 
of the 5-d D^(0) are seen when the 4-d D ov (0) has zero eigenvalues. However, as shown 
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in Sec. [II D| , the pseudo-fermion term has a zero eigenvalue in the infinite fifth dimensional 
extent limit at the zero crossings of H w . We show in Fig. |] the lowest eigenvalue of the her- 
mitian version of the unprojected D^ 5 '(0) and the pseudo-fermion operator D^(l) (recall 
that the pseudo-fermion operator is unaffected by projection). The pseudo-fermion eigen- 
values track the regular eigenvalues up to the crossing, then move away from zero again. 
As L s increases, the pseudo-fermion eigenvalues appear to go exponentially to zero. We 
see then by studying eq. ( JT7\ ) that the pseudo-fermions cancel to some extent in the ratio 
(D^il^D^iO). 

For the resultant unprojected 4-d operator D tov (0) (finite L s ), the eigenvalues are larger 
than the corresponding 5-d ones and the 4-d eigenvalues decrease across a broad region in 
the mass M. These eigenvalues are the actual ones occurring in 4-d observables like the pion 
propagator and the chiral condensate. More to the point, the eigenvalues of the 5-d operator 
are in general not directly physically relevant to 4-d observables. Stated differently, 
the eigenvectors or some piece of the eigenvectors of the D^ 5 '(m) are not eigenvectors of 
Dtovijn) because of the basis changing matrices V and the projection onto the first 5-d 
slice. However, eq. ( |77| ) gives us the connection between the 5-d and 4-d matrices and in 
fact provides a new algorithm to compute the relevant 4-d eigenvalues from only the 5-d 
operators. It has some attractive features for practical use because of its simplicity and 
direct connection, but it is not really efficient for the reason that D^(l) has near zero 
eigenvalues at H w zero crossings and hence has bad condition numbers, at least in the large 
L s limit. 

We proceed to another test of chiral symmetry given by the generalized Gellmann-Oakes- 
Renner relation (called GMOR by many authors). In the form used here and with our 
normalization conventions |J, it states that 

m (b\( l5 D^(m)) 2 \b) = (b\D-\m)\b) V 75 |&> = ±\b) , (80) 

with the physical (subtracted, see eq. (fi5|)) fermion propagator 

D^(m) = - J— [D-l(m) - l] (81) 

1 — TO 

and L s taken to the infinite fifth dimensional limit. Averaging over several (chiral) Gaussian 
random vectors \b), eq. (|30"D becomes a stochastic estimate for mx-K = (V^), which is the 
familiar Gellmann-Oakes-Renner relation. Eq. fl8"0p holds configuration by configuration for 
any chiral state \b). For finite fifth dimensional extent L s , the relation is broken and a useful 
measure is the ratio 



(b \(l 5 Dr o l(m)) 2 \ b) 
(b\D£(m)\b) 



R = m ^, 5 ~T • (82) 



This ratio has been studied for domain wall fermions where in quenched SU(3) configu- 
rations at /3 = 5.85 and 5.7 significant deviations of the ratio from one are seen for fermion 
masses on the order m = 0.01 and smaller. We argue that the larger violations seen at 
P = 5.7 are due to the larger ensemble averaged density of zero eigenvalues of H w (—M) at 



(3 = 5.7 compared to (3 — 5.85 [13 
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We show in Fig. |5] a plot of the ratio R for the domain wall operator as a function of the 
fermion mass at a fixed M = 0.5 and = 1 with the instanton background used in Fig. [l|. 
The mass M = 0.5 is just before the zero crossing of H w (—M). Without projection, the 
deviation is large even for L s = 32. With 5 eigenvalues projected, the deviation from unity 
is negligible for L s = 16. For L s = 8, the violation is large with 5 eigenvectors projected 
out and is not noticeably changed with 20 eigenvectors projected out. This is explained by 
noticing in the right panel of Fig. |3| that for L s = 8 the fifth eigenvalue shows a deviation of 
1 — |^4 (A5) I ~ 0.01. What is not shown is that this value changes little at the 20th eigenvalue. 
What we are seeing is that by the 20th eigenvalue we are entering into a relatively dense band 
of eigenvalues of and they change slowly with increasing eigenvalue. Therefore, little is 
gained by projecting out more eigenvectors. By going to L s = 16, 1 — ( A5) | ~5x 10~ 5 is 
much reduced compared to L s = 8, and hence all the volume modes contributing to R have 
deviation less (or much less) than the worst case resulting in a small overall deviation of 
R from unity. Therefore, projection is quite effective at restoring chiral symmetry, and we 



see that the GMOR test given by eq. (|82|) is a quite sensitive measure of chiral symmetry 
breaking since it detected a deviation in the e approximation on order of 1%. 

We finish up the discussion of the smooth field case by considering the case of M > 2 
(the doubler region). We continue the discussion initiated in Ref. |T7| on the difference of 



the overlap Dirac operator D ov (m\H w ) and the 4-d version of the standard domain wall 
operator D ov (m; H T ). 

We consider the free propagator with momentum near one of the corners of the Brillouin 
zone, i.e., with p^ = sin(p M ) <C 1 for all /i. Set M{p) = M — B{p) = M — ^ 2 sin 2 (p M /2) = 
M— 2n+0(p /J 2 ), where n is the number of momentum components near 71". Then D w (—M) = 

1 ^\ IfiPfi ~ M{jp) and we find for the inverse free overlap propagator, up to terms of higher 
order in and m, 

(2M(p)x ~^l^ +2 ZrS forM(p)>0 
l^ir forM(p)<0. 

This corresponds to a light free fermion propagator only for M(p) m M — 2n > 0. For 

2 < M < 4, for example, the origin and the corners of the Brillouin zone with one momentum 
component close to ir give light fermions, while for < M < 2 only the origin gives a light 
fermion. 

For free domain wall fermions, recalling the form of H? from eq. (EO) , 



Ht{ . M) = 75fl „ ( _ M) _L_ = 75fl „ ( _ M) _l_ . (84) 



we find 

D T (-M) = l5 H T (-M) = C(p) 



Y^1,P,-\m{ P ) (2-a 5 M(p)) 



i 



(85) 



with C(p) positive, and hence unimportant in e(Hj<)- We see that for free domain wall 
fermions |M(p) (2 — a^M{p)) plays the role of M(p) for overlap fermions. Hence, from 
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eq. (^3|) we see that we need M(p) (2 — a^M{p)) > to have a light domain wall fermion. 
These are exactly the conditions for a normalizable domain wall zero mode jLTj], —1 < 
1 — a^,M(p) < 1. The first inequality here is an additional condition for domain wall 
fermions, not present for overlap fermions. For 2 < M < 4 it excludes the origin of the 
Brillouin zone from giving a light fermion, as it did for overlap fermions. 

Now consider the interacting case. We know that zero crossings of H w (—M) = also 
correspond to zero crossings of the domain wall Ht{—M) for all 05. However, eq. (|S4J) 
suggests that H T (—M) can have poles. Zero eigenvalues of H w (—M) at some M are also zero 
eigenvalues of D w (— M ). From (IS1) we see that to each such zero eigenvalue is associated a 
pole of Hj-(—M) at M = M + 2/0,5. 111 other words, all zeros below M = 2 are replicated 
as poles at M + 2/a^. 

For the single instanton background shown in Fig. [I], there is one downward crossing in 
the Wilson spectral flow at about M = 0.55. After M = 0.55, the domain wall and overlap 
Dirac operators have index —1 and hence a single zero mode. For this one downward 
crossing, there are 4 corresponding upward crossings at M = 2.17 and M = 2.30. After this 
last crossing, the index of the Dirac operator is +3 and hence there are 3 zero modes. There 
are no other zeros of H w (—M) and hence Ht(—M) until M ~ 4. The largest eigenvalue of 
Ht(—M) has a behavior as follows: Starting at M — 0, there is a near pairing of positive and 
negative eigenvalues of the largest eigenvalue A max ; however, a negative eigenvalue is slightly 
largest in magnitude. As M increases, so does |A max |. At M — 1.249 A max reaches —1, giving 
a zero eigenvalue for T~ l . This corresponds to B in the transfer matrix T, eqs. (|i~7j ) and 
(UP, having a zero eigenvalue. In the free field case the domain wall mass at which T~ l has 
a zero eigenvalue (namely M = 1) is identified as the mass where the domain walls (each end 
of the fifth dimension) have the least coupling for fixed L s , i.e. this is the optimal domain 
wall mass. As M increases to the pole position of M = 2.55, the eigenvalue — A max diverges. 
The pole is the mechanism by which the topology changes for as M increases beyond the 
pole mass the largest eigenvalue flips sign and decreases from positive infinity. There is a 
net increase in the number of positive states of Ht{—M) occurring in e(Hx(—M)) for the 
4-d domain wall Dirac operator D ov (0; Hj<(—M)) and the index becomes +4 with four zero 
modes. 

In general then for the standard domain wall action in the smooth field case with a$ = 1, 
for every crossing of H w (—Mq) for Mq < 2, there are four opposite crossings around M > 2 
and a pole of Ht(—M) at M = 2 + Mq. If at some mass M < 2 there is an index Q(M) 
then Q(2 + M) = -AQ{M). For the overlap Dirac operator D ov (0; H w (—M)), however, 
Q(2 + M) = -3Q(M). 

Varying the fifth dimension lattice spacing 05 gives some freedom in improving the chiral 
symmetry properties for the induced 4-d Dirac operators at finite L s given in eq. fl78|) . If at 
some positive A* < 1 we have 1 — £l 3 /2(<25A*) < 5 where 5 is some prescribed accuracy, then 
the range of accuracy is |A| in [A*, l/(a§A*)]. For the overlap Dirac operator with H w (—M), 
the flexibility of rescaling by a gauge field dependent a 5 coupled with projection helped 
insure that all unprojected eigenvalues of H w (—M) where within some prescribed range of 
accuracy of \eL s /2(cL 5 \i)\ approximating 1. 

For the standard domain wall operator, Ht{—M) depends implicitly on as, and choosing 
different values for 05 leads to different cut-off effects, just as choosing different values for M 
does. Therefore, both M and a 5 need to be kept fixed in a given simulation. We saw in Fig. |3| 
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that even for moderate L s the largest eigenvalue A max of Ht(—M) was usually well within 
the good approximation region of e(a 5 A). Thus we would like to use an a 5 > 1 to fully use 
the region of good approximation to e. However, this can mix in doubler contributions, since 
for each crossing at M < 2 there is the doubler pole of H T (—M) at M = 2/a 5 + M , and as 
as increases the pole moves down to lower domain wall mass M. However, as can be seen 
in the right panel of Fig. |3] the deviation for the largest eigenvalue of Ht becomes sizeable 
well before the pole at M = + Mq is reached. This region of large deviation moves 
down with increasing a$ and begins to conflict with our goal of increasing to fully utilize 
the region of good approximation to e. Since one needs to choose M sufficiently larger than 
the M c around which the physically important crossings, due to large instantons, occur, the 
choice of as is quite limited. This is a fundamental limitation of the standard domain wall 
operator. 

B. Quenched gauge field 

We now move on to the interesting case of a thermalized quenched SU(3) gauge field 
configuration. Fig. |] shows the spectral flow of the hermitian domain wall and Wilson 
operators for a quenched SU(3) (3 = 5.85, 8 3 x 32 configuration. Near M = 1.0, the Wilson 
flow has two crossings down and two crossings up that result in a zero index for M > 1.017. 
There is a fifth crossing at about Mq = 1.86 resulting in an index of —1 for M greater 
than this crossing mass. In the region around M = 1.0 where there is a non-zero index, 
the domain wall eigenvalues do not go close to zero. Well into the region M > 1.017 where 
the index is zero there are fairly small domain wall eigenvalues, but no dramatic drop of 
the lowest eigenvalue is observed for M > Mq. The overlap Dirac operator, on the other 
hand, has zero modes whenever the index is non-zero. Fig. shows the spectral flow of 
the projected domain wall operator on the same configuration as in the Fig. |^, together 
with the Wilson flow and the low-lying eigenvalues of 2Ht- Near M = 1.0 and after the 
Wilson crossing around 1.85 there are small (zero) eigenvalues of the projected domain wall 
operator. 

In the left panel of Fig. |8| is shown the spectral flow of the lowest eigenvalue of the 
standard domain wall operator for L s = 10 and 30 and the Borigi variant for L s = 10 all for 
05 = 1. We see in Fig. |8| that a small (near) zero eigenvalue is obtained for L s = 30 with 
the standard domain wall D^ w (0). For L s = 10, a drop is seen in the lowest eigenvalue for 
M > Mq, down to the level of the L s = 30 eigenvalue. For the Borigi variant Dp W ,(0), a 
relatively small eigenvalue is seen even for L s = 10. 

The explanation for the small eigenvalues of the Borigi variant even at L s = 10 can be seen 
in the right panel of Fig. || which shows the accuracy of the |ez, s /2(Aj)| approximating 1. The 
lowest eigenvalues of Ht(—M) are about half those of H w (—M). Hence for the domain wall 
fermion action with L s = 10 and M = 1.65 we find 1 — |e^5 ( A20) I ~ 0.25 while for the Borigi 
variant using H w (—M) the deviation is about 0.3%. The largest eigenvalue of Ht(— 1.65) is 
about 1.13 hence the deviation is below the range in the graph even for L s = 10. For the 
Borigi variant, the largest eigenvalue of H w (— 1.65) is 5.82 and there is about 6% deviation. 
While the condition number is smaller for Ht{— 1.65) than H w {— 1.65), the eigenvalues are 
placed more optimally around one for H w which results in overall smaller deviations of 
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\ £ L s /2{^i)\ from 1. For this configuration, there is not much need for rescaling the Wilson 
H w eigenvalues by 05 since 05 = 1 is close to optimal. For the standard domain wall action, 
a larger 05 would lower the deviations of |£l s /2(g5A;)| from 1 for this configuration. However, 
as explained in the previous sub-section, choosing an a 5 > 1 might not be desirable, since 
it could lead to large deviations (due to the proximity of poles) for other configurations in 
a simulation. 

Fig. ^ shows that small eigenvalues of the domain wall fermion actions D^ w and Dp W , 
can be easily achieved using projection and the cost of the projection is justified. To achieve 
the same small eigenvalue without projection would require L s ^> 30. 

In Fig. |9| is shown the ratio R defined from the generalized Gellmann-Oakes-Renner 
relation (GMOR) in ([82]) using M = 1.65, as = 1 and the same configuration as in Fig. || 
For the standard domain wall action, not much is gained from projection for L s = 10 and 
violations of the relation are large. For L s = 30, projection significantly improves chiral 
symmetry resulting in small violations of the relation while there are large deviations for 
even this large L s without projection. The Borigi variant also shows similar improvement 
with increasing L s . 

The condition number of a fermionic operator is a useful measure of its convergence 
properties in solving linear systems of equations, and can be used to construct a bound on 
the number of iterations for convergence in some iterative methods. The condition number 
defined here is A max /A min where Aj > is defined from D^Dvi = Xfvi. The condition numbers 
of the various fermion operators are shown in Fig. pi]. The parameters are chosen the same 
as for the GMOR test in Fig. || The condition number of the overlap Dirac operator 

D ov (m; H w ) is much smaller than the 5-d methods. The condition number of the standard 

(5) 

domain wall operator is larger for L s = 30 compared to L s = 10 while the Borigi 

(5) 

variant is roughly the same for both L s values, but higher than that of D DW . Projection 
also does not affect the condition numbers much basically because for this configuration at 
mass M = 1.65, there should be no zero modes and the smallest eigenvalue of D^> is not 
strongly affected by projection as can be seen in Fig. [| The reason the 4-d operator has 
much smaller condition number is due to the largest eigenvalue which for D ov (m) is very 
close to 1, while for it is roughly 10. 



Performance tests and comparisons of the various fermion actions are shown in Fig. ITT 



The left panel shows the iteration count for the benchmark Conjugate-Gradient for Normal 
Equations (CGNE) algorithm to achieve a residual accuracy of 10~ 5 normalized by the 
source norm for the linear system D^Dcf) = D^b. The configuration parameters are those of 
the GMOR test in Fig. []. The average and the standard error of the data (not the mean) 
is derived from the 12 source color and spin inversions necessary to compute a fermion 
propagator. For the 4-d case (exact chiral symmetry not necessary), we note that a single 
inversion could be used for all the fermion masses m used in Fig. [TTJ with the convergence 
governed by the time required for the smallest fermion mass. The 4-d overlap Dirac operator 
D ov (m) has a low iteration count requiring about 320 CG iterations for convergence at m = 0. 
Without preconditioning, the standard domain wall operator D^ w (m) requires roughly 1000 
iterations for convergence at L s = 10 and there is a strong L s dependence as expected from 
the condition number. Preconditioning reduces the condition numbers by about a factor 
of three. The Borigi variant requires correspondingly more iterations for convergence. All 
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these results are consistent with a linear scaling in the condition number. 

It has been pointed out 0,0] that both CGNE and Conjugate- Residual are optimal 
algorithms for the overlap-Dirac operator with exact chiral symmetry. In principle the 
Conjugate- Residual method would be more desirable since only one D ov {m) application is 
needed per iteration. For CGNE one expects two applications of D ov (m), but this can be 
rewritten as only one application of e(H w ) on a vector fl2fl , so the work is the same per 
iteration for both methods. However, it was found that Conjugate- Residual had clearly 
worse convergence in practice compared to CGNE and was not used. 

The right panel of Fig. [TT] shows the actual cost for convergence measured in total number 
of D w applications for each linear system solution. This includes the "inner" CG portion 
for the 4-d overlap Dirac operator. For each "outer" CG iteration, a single application of 
t{H w ) is needed on a vector ||. This inner CG iteration in turn involves multiplications of 
a vector with H^. For the standard domain wall operator, there are 2L S applications of D w 
needed per iteration. For the Borigi variant, there are 3L S applications of D w needed as can 
be seen from the form of Z)^^,(m) and (D^^^m))^ in eqs. (|10|) and (|58|). 

The cost of the methods can be roughly characterized as follows: the Borigi variant 
requires about twice as many D w applications for a given L s compared to the standard 
domain wall operator with and without projection. This can be expected from the slightly 
larger condition number (mainly from the largest eigenvalue) and the 50% overhead in D w 
applications. Preconditioning saves the domain wall operator roughly a factor of two to 
three for both L s values with and without projection. This decrease is mainly due to the 
reduction of the largest eigenvalue. The L s = 30 calculations require about three to four 
times as many D w applications as the L s = 10 calculations with and without projection. 
However, L s = 30 and projection is needed to pass the GMOR test successfully. 

The inversion of (1 — A ee P_) in eq. (|7^) requires only about three iterations of a Min- 
imal Residual algorithm to reach an accuracy of 10~ 7 . This multiplicative overhead in the 
number of "outer" CG iterations, but not proportional to L s , is included in the total cost of 
the preconditioned domain wall operator and results in about an 18% overhead for L s = 10 
and a 5% overhead for L s = 30 — there is little L s dependence in the "inner" CG. However, 
for these parameters and configuration (topological index is zero) the projected precondi- 
tioned domain wall requires fewer CGNE iterations for convergence, and is less costly in D w 
applications for L s = 30 than the unprojected version. 

For L s = 30 with preconditioning and projection, the standard domain wall operator 
is about three times less costly in D w applications than the 4-d overlap Dirac operator. 
However, given that with the 4-d overlap Dirac operator all the masses for this spectroscopic 
calculation can be computed simultaneously for a given color and spin source, the cost benefit 
of the two approaches is comparable. 

IV. DISCUSSION 

The five dimensional domain wall fermion action provides a means whereby an effective 
chiral theory may be obtained. Upon integrating out extra modes (which may be interpreted 
as flavors) from this extra fifth dimension, a chiral theory is obtained in four dimensions. 
For finite fifth dimensional extent L s , the form of this action in four dimensions is 
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D tov (m;H T ) = - 



1 + m+ (1 - m)%e Ls/2 (a 5 H T ) (86) 



with the auxiliary Hamiltonian H T . When L s — > oo, eL g / 2 (a 5 H T ) — > e(H T ) for a 5 > 
and there is a chiral symmetry in the action. For finite L s , chiral symmetry breaking is 
induced because 6L a / 2 (a 5 H T ) deviates from one. These deviations occur whenever the range 
of eigenvalues of H T are outside the range of approximation of Sl s / 2 (x) to e(x). Deviations 
can occur for both large and small eigenvalues of Ht- 

Section || shows how the domain wall operator and a variant can be straightforwardly 
modified to allow arbitrarily precise chiral symmetry at finite values of L s . The domain 
wall operator as shown in eq. (|T0|) has only two new entries corresponding to projection of 
eigenvalues of Ht or H w . When the \ basis is chosen for the fermions fields in the action, 
eq. (|2T|), the projection term occurs only in the origin along the fifth dimensional line (the 
flavor index). With this special position, the projection term is unmodified by the successive 
integrations of the extra flavor fermion fields. After extraction of the pseudo-fermion term, 
the expected form of the 4-d fermion action is obtained with the projection term acting as 
(part) of the spectral decomposition of e{a^H). In principle then one can see that once in 
the x basis the projection term could account for the entire spectrum of H and the fifth 
dimensional extent could be set to one! 

All the normal relations involving the fermion propagator go through with projection, 



and we have the relation connecting the 5-d and 4-d operators in eq. ([51]) and the propagator 
in eq. (^) valid with and without projection. 

The main purpose of this paper is to show how projection can be used to achieve exact 
chiral symmetry at finite fifth dimensional extent. The method is not limited to quenched 



calculations and methods were outlined in Section |II F| for use in dynamical fermion calcula- 
tions. Even-odd preconditioning is practical and for moderate L s is a small overhead. Tests 
of chiral symmetry were made for a smooth gauged field background composed of a single 
instanton, and for a quenched SU(3) gauge field background. Small (zero) eigenvalues of the 
Dirac operator can be achieved given careful control of the range of the approximation of 
the eigenvalues of the auxiliary Hamiltonian and choice of L s . A stringent simple test is the 
generalized Gellmann-Oakes-Renner relation (GMOR) from eq. fl8~0|) and projection appears 
essential to pass this test. 

However, for practical simulations the requirement of no induced chiral symmetry break- 
ing may not be necessary or easy to achieve. In this case, the amount of projection and the 
fifth dimensional length can be tuned. This may well be necessary anyway since as shown 
for quenched backgrounds, there is evidence for a non-zero density of zero eigenvalues 
p(0; H) of the auxiliary Hamiltonian H w (and hence Ht) for a large range of (3. If one wants 
to go to very large lattice volume, the number of eigenvalues for projection to achieve some 
given accuracy for EL a / 2 {a§H) grows like the 4-d lattice volume. For example, this cost is 
an additional overhead for a spectroscopy calculation which also grows like the 4-d lattice 
volume. One way to lessen the overhead is to use a weaker coupling since p(0; H) decreases 
very rapidly with the coupling which of course may require a larger lattice volume to hold 
the physical volume fixed. To keep calculation costs down, L s can be chosen suitably small 
and measures of induced chiral symmetry breaking like PCAC [0] can be used. No effort 
was made to quantify this for given parameters like however, some projection is expected 
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to reduce the induced chiral symmetry breaking from finite L s . 

For a given calculation, one would like to know what ultimately is the most efficient 
method to use. Namely, is the 5-d or the direct 4-d method the most efficient? From the 
present work, when one enforces chiral symmetry there is not really a huge difference in the 
cost of any of the methods tested. Quite likely then the answer depends on what one wants to 
do and on ones taste. For eigenvalue calculations of spectral quantities, the 4-d eigenvalues 
are needed. For spectroscopy or dynamical fermion calculations, either 5-d or 4-d methods 
can be used. For fixed cutoff dependence, the 5-d domain wall operator is preferable to 
the 4-d methods using e{a^,Hj) - the nontrivial cost of applying Ht within an inner CG 
iteration inside of another inner CG iteration is prohibitive. For 4-d eigenvalue calculations 
the application of D ov (m; Ht) from D^ w (m) in eq. (|5l"D is probably most efficient. For 
the overlap Dirac case using D ov (m; H w ), the Borigi variant of D^' w ,(jn) does not appear 
competitive for spectroscopy calculations. A multi-grid implementation in the fifth direction 



might help [I4"|, but has not been tested here. Other 5-d methods should be tried and 
tested [^3] . 



For a direct comparison of the standard domain wall action D^ w (m) and the overlap 
operator D ov (m), one first must answer just how important is chiral symmetry? Projection 
is essential for realizing good chiral symmetry properties for small quark masses and mod- 
erate L s . Clearly computing eigenvectors using H w is much more efficient than computing 
eigenvectors of Ht- If one enforces chiral symmetry the un-preconditioned domain wall op- 
erators is about the same cost in the number of D w applications as the overlap operator. 
Preconditioning saves about a factor of three in cost, however the methods are comparable 
in cost for a spectroscopy calculation when multiple fermion masses are needed. A benefit 
of the 5-d method is that it gives a simple handle L s in which to reduce the overall cost 
while incurring some acceptably small amount of induced chiral symmetry breaking. 
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FIG. 1. Spectral flow of the hermitian DWF operator for various L s and of the hermitian Wilson 
operator around its zero crossing on a single 8 4 instanton background. Also shown is the flow of 
the overlap-Dirac operator which has zero eigenvalues after the crossing, setting in discontinuously. 
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FIG. 2. Spectral flow of the projected DWF operator, projecting the five lowest eigenvectors 
of Ht, for various L s on the same instanton background as in Fig. [I]. With projection, the lowest 
DWF eigenvalue drops to (near) zero discontinuously after the Wilson crossing. Also shown are 
the flows of H w and 2 * Ht- The eigenvalues of 2 * Ht track closely those of H w . 
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FIG. 3. Results for DWF on the same 8 instanton background as in in Fig. |. The left 
panel shows a logarithmic plot of how the lowest (near) zero eigenvalue (with extent L s and iV 
eigenvectors projected) slowly rises with decreasing accuracy of the epsilon approximation as shown 
on the right. The right panel shows the quality of the approximation of |ei, s /2(\)| to unity for 
eigenvalues Aj of Ht- The Ai deviations shows the quality for the standard action, Atv is for 
the fifth eigenvalue, and X max is for the highest eigenvalue. The smallest projected domain wall 
eigenvalue rises as the deviation for X max increases. 



29 



0.4 



0.3 



< 0.2 



L s = 8: o m 

= 16: 



L 



L 



= 0. 



1X1 : 



u. 




□ m 

□ m 

□ m 



L s = 24: o m 
L s = 32: o m=0, □ ni ; 
s = 48: o m=0, □ m : 
- Wilson 



= 1 
1 
1 
1 
1 



0.1 




□ 
o 



0.4 



0.6 



0.8 



FIG. 4. Spectral flow of the lowest pseudo-fermion eigenvalue (m = 1) and lowest domain wall 
eigenvalue (m = 0, with no projection) for various L s on a single 8 instanton background. The 
minimum of the lowest pseudo-fermion eigenvalues is at the Wilson crossing and decreases with 
increasing L s . 
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FIG. 5. The ratio given by the Gellmann-Oakes-Renner relation (GMOR) for M = 0.5 (just 
before the crossing) as a function of L s and projection on a single 8 4 instanton background. N 
denotes the number of projected eigenvectors. Projection for L s = 16 helps to satisfy GMOR since 
for all unprojected eigenvalues of Ht, e(-Hr) is well approximated. For L s = 8, the deviation for 
the largest projected eigenvalue is still large as seen in Fig. |3[ 
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FIG. 6. Spectral flow of the hermitian (unprojected) DWF operator for various L s , and of the 
hermitian Wilson and overlap-Dirac operators on a quenched SU(3), (3 = 5.85, 8 3 x 32 configuration 
around the zero crossing of the Wilson flow. The overlap operator correctly finds zero eigenvalues 
around M = 1 while the domain wall operator misses them. 
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FIG. 7. Spectral flow of the projected DWF operator, projecting the 20 lowest eigenvectors of 
Ht, for various L s on the quenched SU(3) configuration of Fig. |6[ With projection, the lowest 
domain wall eigenvalue drops to zero (within precision used here) discontinuously after the crossing 
around 1.85. The zero eigenvalues around M = 1 are also correctly identified. Also shown are the 
flows of H w and 2 * Ht- The eigenvalues of 2 * Ht track closely those of H w . 
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FIG. 8. Results for the quenched SU(3) configuration of Fig. ||. The presentation is similar 
to Fig. ||. DWF and OVDWF are the domain wall operator and Borici's variant, resp., with N 
eigenvectors projected. The left panel shows how the lowest (near) zero eigenvalue slowly rises with 
decreasing accuracy of the epsilon approximation as shown on the right. The right panel shows 
the quality of the approximation of \^L s /2{^i)\ f° r eigenvalues A, of Ht- Here N = 20 eigenvectors 
are projected. For domain wall, the 20-th eigenvalue is badly approximated for L s = 8 but slightly 
better for the Borigi variant. For L s = 30, the deviation of A max is below the range shown. 
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FIG. 9. The ratio given by the Gellmann-Oakes-Renner relation (GMOR, c.f. (|80|) ) for 
M = 1.65 (just before the crossing) as a function of L s and projection on the same quenched 
SU(3) configuration as in Fig. |(| DWF and OVDWF are the domain wall operator and Borigi's 
variant, resp. For L s = 30 with projection, the GMOR relation is well satisfied but some discrep- 
ancy is seen from deviations of Ht eigenvalues Aj in Fig. ||. 
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FIG. 10. Condition numbers as a function of the mass for various actions on the same quenched 
SU(3) configuration and parameters as in Fig. |9[ The overlap operator has the lowest condition 
number. 
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FIG. 11. The left panel shows the number of (outer) CGNE applications to reach convergence 
of 10 -5 . The right panel shows the total number of D w applications for convergence (including 
inner and outer CG applications) on the same quenched SU(3) configuration and parameters as 
in Fig. |9|. The overlap operator needs the fewest CG iterations. For L s = 30 with projection, the 
standard (un-preconditioned) domain wall and overlap operator have comparable D w application 
counts. Preconditioning saves the domain wall operator about a factor of two to three in cost. 
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